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Abstract 

We study the dynamics of systems on networks from a linear algebraic perspective. 
The control theoretic concept of controllability describes the set of states that can 
be reached for these systems. Under appropriate conditions, there is a connection 
between the quantum (Lie theoretic) property of controllability and the linear systems 
(Kalman) controllability condition. We investigate how the graph theoretic concept of 
a zero forcing set impacts the controllability property In particular, we prove that if 
a set of vertices is a zero forcing set, the associated dynamical system is controllable. 
The results open up the possibility of further exploiting the analogy between networks, 
linear control systems theory, and quantum systems Lie algebraic theory. This study 
is motivated by several quantum systems currently under study, including continuous 
quantum walks modeling transport phenomena. Additionally, it proposes zero forcing 
as a new notion in the analysis of complex networks. 

1 Introduction 

This paper deals with several concepts from different fields such as linear algebra, 
graph theory and quantum and classical (linear) control theory. In the context of 
dynamics and control of systems on networks, it establishes a connection between a 
notion of graph theory (zero forcing) and concepts in control theory (quantum and 
classical controllability). We review these different concepts before we introduce the 
technical content of the paper and give physical motivation for our study. 
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1.1 Background 



For a dynamical system with a control input, the property of controllability describes to 
what extent one can go from one state to another with the evolution corresponding to 
an appropriate choice of the control. If all the possible state transfers can be obtained 
within a natural set (the phase space), then the system is said to be controllable. 

For several classes of systems, controllability has been described in detail and 
controllability tests are known. In particular, for linear systems 

s 

x = Ax + ^b i% , (1) 

i=i 

A 6 R" x ", bj e K n , j = 1,2,..., s, where both the state x 6 M" and the control 
functions Uj = uj(t) enter the right hand side linearly, several equivalent conditions 
of controllability are known. The classical Kalman controllability condition (see, e.g., 
[TS]) says the system ([1]) is controllable if and only if the n x (ns) matrix 

W(A, B) := [bi, Abi, . . . , A n ~ 1 bi, ....Ik., lb A n_1 b s ], 

has full rank n, where B :~ [bi \>2 ■ ■ -b s ]. In this case, for any prescribed state 
transfer Xo — s- xi(s W 1 ) and interval [0,T], there exists a control u(t) = [ui, . . . ,u s ] T 
such that the corresponding solution x(t) of ([T]) satisfies x(0) = xo and x(T) = 
xi. For quantum mechanical systems which are closed (i.e., not interacting with the 
environment) and finite dimensional, one considers the Schrodinger equation 

ij t \^=H(uM, (2) 

where \ip) G C" is the quantum state and the Hamiltonian matrix H = H(u) is 
Hcrmitian and depends on a control u = u(t) which in some cases can be assumed 
to be a switch between different Hamiltonians. If is a system linear in the state 
|V>), the solution of © is \ijj(t)) = X(t)\ip(0)) where X = X(t) is the solution of the 
Schrodinger matrix equation 

iX = H(u)X (3) 

with initial condition equal to the n x n identity matrix /„. Since H = H(u) is 
Hermitian for every value of u and therefore —iH is skew-Hermitian, the solution 
of ([3]) is forced to be unitary at every time t. In this context, the system is called 
completely controllable if for any unitary matrix Xf in S*?7(n|3 there exists a control 
function u = u(t) and an interval [0, T] such that the corresponding solution X = X(t) 
of © satisfies X(0) = I n and X{T) = X f . 

At the beginning of the development of the theory of quantum control, it was 
realized (see e.g., [TT]) that system ((3]) has a structure familiar in geometric control 
theory |13] and therefore controllability conditions developed there can be directly 
applied. In particular, the Lie algebra rank condition [14] says that a necessary and 
sufficient condition for complete controllability of system is that the Lie algebra 
generated by the matrices {iH(u)} (as u varies in the set of admissible values for 
the control) is su(n) or w(n,)H This has given rise to a comprehensive approach to 

1 Following standard notation, SU (n) is the special unitary group, i.e., the matrix group of n x n unitary 
matrices having determinant 1. 

2 Following standard notation, u{n) is the Lie algebra of n x n skew-Hermitian matrices and su(n) is 
the Lie algebra ofmxm skew-Hermitian matrices with zero trace. 
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quantum control based on the application of techniques of Lie algebras and Lie group 
theory [7]. 

In recent years there has been considerable interest in the study of control systems, 
both classical and quantum, which are naturally modeled on networks. Often one tries 
to relate the controllability of these systems to topological or graph theoretic proper- 
ties of the network. For quantum systems, the nodes of the network may represent 
energy levels or particles which are interacting with each other. For these systems, the 
application of the Lie algebra rank condition to determine controllability can become 
cumbersome and subject to errors when the dimension of the system becomes large. 
It is preferable to have criteria based on graph theoretic properties of the network not 
only because they are typically checked more efficiently but also because they give 
more insight in the dynamics of the system. Work in this direction has been done in 
[2], 0, [IB]- In this context, a relevant property of a graph G and a subset S of its 
vertices is the capability of this set to 'infect' all the vertices of the graph, as explained 
in the next paragraph. 

Every graph discussed is simple (no loops or multiple edges), undirected, and has 
a finite nonempty vertex set. Consider a graph G and color each of its vertices black 
or white. A vertex v is said to infect, or force a vertex w if v is black, w is white, w is a 
neighbor of v, and w is the only white neighbor of v. In the case where infection of w 
has occurred, we change the color of w to black and continue the iterative procedure. 
The set S is called a zero forcing set if this procedure, starting from a graph where only 
the vertices in 5" are black, leads to a graph where all vertices are black. An example 
of a zero forcing (infection) process is shown in Figure [TJ indicated by arrows; the set 
of black vertices is a zero forcing set. 
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Figure 1: A zero forcing set and the process by which it can infect all vertices. 



For a real symmetric n x n matrix A = [dkj], the graph of A, denoted G(A), is 
the graph with vertices {1, . . . , n} and edges {kj : a^j ^ and k ^ j}. Observe that 
G = Q(Ag) = G(Lg), where Ac and Lq = Dg — Ac denote the adjacency matrix 
of G and the Laplacian matrix of G, respectively (here Dg is the diagonal matrix 
of degrees) . Zero forcing has been studied in detail in the context of linear algebra. 
This is because the size of the minimum zero forcing set of a given graph G, which 
is called the zero forcing number Z(G), is an upper bound to the maximum nullity 
(or maximum co-rank) over any field of G [3]; the maximum nullity is taken over all 
symmetric matrices A such that G(A) = G (see [8] for background on the problem of 
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determining maximum nullity). 

Zero forcing appears then to be a valuable concept in the study of graph-theoretic 
properties that are captured by generalized adjacency matrices. Indeed, there are im- 
portant classical parameters introduced with this purpose, e.g., the Colin de Verdierc 
number, the Haemcrs bound, etc. It has to be remarked that questions about the max- 
imum nullity of a graph are generally difficult problems and the zero forcing number 
does not constitute an exception: it was shown in [T] that there is no poly-logarithmic 
approximation algorithm for the zero forcing number. 

1.2 Contribution of the paper and physical motivation 

In this paper, we consider the dynamics of a system defined on a network and relate 
the above notions and criteria of controllability with the graph theoretic concept of 
zero forcing. Abstractly, we consider a graph G and a subset S = {ji, ■ ■ ■ , j s } of 
its vertices V(G) = {1, ...,n}. The dynamics are that of a quantum system <|3j> 
where the Hamiltonian is allowed to take the values {A, e^e^, . . . ,e Js eJ s }. Here A 
is the adjacency matrix Aq of G, Laplacian matrix Lq of G, or more generally a 
real symmetric matrix such that G(A) = G with all nonzero off-diagonal entries of 
A having the same sign (which is the typical situation in transport models). The 
vectors {e^ , . . . , e^} are the characteristic vectors^ of the vertices in S. In this way, 
we can associate a linear system ([1} with A and hi = , . . . , b s = ej 3 . The main 
result of the present paper says that controllability in the quantum sense, expressed 
by the Lie algebra rank condition, and controllability in the sense of linear systems, 
expressed by the Kalman rank condition, are equivalent conditions. Moreover, if 
the set S (corresponding to e^, . . . ,Bj a ) is a zero forcing set, then these equivalent 
controllability conditions are true (the converse is false). The first of these results is 
along the same lines as the main result of |10j which considers the case of quantum 
dynamics switching between the Hamiltonian A and zz T , where z = Yljes e o> ano - 
establishes the connection between controllability (quantum and linear) . As mentioned 
above, these characterizations avoid lengthy calculations of the Lie algebra generated 
by a given set of Hamiltonians and replace them with more easily verified graph 
theoretic and linear algebra tests. 

On physical grounds, our motivation for considering a Hamiltonian specified by a 
matrix with the given graph comes from the study of continuous time quantum walks 
which model transport phenomena in many physical and biological systems [6]. A 
recent review is given in [4]. Most of the studies consider this sole Hamiltonian and 
concern the statistical (diffusion) properties of the dynamics. We add here the Hamil- 
tonians ejej where ej is the characteristic vector of a given node of the network and 
study the nature of the states that the resulting dynamics can achieve, in particular 
whether an arbitrary (unitary) state transfer can be achieved between the states of 
the quantum system. The Hamiltonians e^ej model a prescribed energy difference 
between the corresponding node and all the other nodes of the network which are 
assumed to be at the same energy level. Therefore the dynamics is the alternating 
of a diffusion process (modeled by the Hamiltonian A) and a rearrangement of the 
energies of the various states by selecting one of the states as high energy state and 
all the other at the same (lower) energy 

Theoretical research in network theory has focused on a number of discrete time, 
deterministic diffusion processes on graphs. While zero forcing has not been studied 

3 The vector &j has the jth entry equal to one and every other entry equal to zero and is also called the 
jr'th standard basis vector. 
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in this context, there are two directions of research that are closely related: as it 
was already noted in [TJ, the threshold model introduced for studying influence in 
social networks shares with zero forcing certain issues underlying its computational 
complexity |16j ; the model of complex networks controllability recently proposed in 
|17j also makes a natural use of the Kalman rank condition and it singles out certain 
combinatorial properties to determine when the condition is satisfied. Determining 
whether zero forcing has a place in the metrology of complex networks is a point worth 
further interest. 

The paper is organized as follows. In Section [5] we introduce notation and give 
background and basic results concerning Lie algebras that will be used in the following 
sections. The connection between quantum (Lie algebraic) controllability and the 
Kalman criterion for linear systems is established in Section [3] There we also prove 
the converse of the main result of [10] . The relation with the zero forcing property is 
established in Section |4j and Section [5] contains concluding remarks. 

2 Lie algebra terminology and preliminary results 

Standard material on Lie algebras can be found in [12] . For Ax,...,A k G C™ xn , 
(Ax, ■ ■ ■ , j 4.fe)[.,.] denotes the real Lie algebra generated by Ax, ...,A k under addition, 
real scalar multiplication, and the commutator operation. Let H n (M) denote the real 
vector space of symmetric matrices. For A G H n (M.), the notation A = [akj] means for 
k < j the (k,j) and (J, k) entries of A are both a k j. Observe that A = [akj] G "H n (R) 
can be expressed as 

n 

A = ^a kk e k e k T + ^a kj (e k ej T + eje k T ). 

k=l k<j 

The following proposition is well known (a proof appears in |10j). It provides a link 
between an appropriate Lie algebra of real matrices and the Lie algebra rank condition 
of quantum control theory, thereby allowing us to work with real matrices only. Recall 
that the Lie algebra consisting of all real n x n matrices is denoted by gl(n, R), sl(n, R) 
denotes the Lie algebra of real n x n matrices with zero trace, u(n) denotes the Lie 
algebra of all skcw-Hcrmitian (complex) n x n matrices, and su(n) denotes the Lie 
algebra of all skcw-Hcrmitian (complex) n x n matrices with zero trace. All these Lie 
algebras are considered as vector spaces over the field of real numbers. 

Proposition 2.1. For Ax,...,A k G H n (R), 

(A 1 ,...,A k ) [ . > . ] =gl(n,R) (iA u . . . ,iA k ) W] = u(n). 

The next lemma is used in the proof of Theorem 13.71 in the next section. 

Lemma 2.2. Let A,Bx,...,B s G %„(R), with s > 1. Define C := (A,Bx, . . - ,-B s )[.,.] 
and let C denote the smallest ideal of C that contains B^, i = 1, . . . , s. If C = gl(n,M) 
and tx B k ^ for some B k , then C = gl(n,M). 

Proof. For n = 1 the result is clear, so assume n > 2, C = gl(n,M), and txB k ^ 
for some B k . Observe that C := (A, Bx, . . ■ , -B«)r.,.] is spanned by A and C. Since 
C = gl(n,M), we have 

\gl(n,M),gl(n,R)] = [span(A) + t, span(A) + £} C £. 
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It is known that [gl(n, R), gl(n, R)] = sZ(n,R), because [gl(n,M),gl(n,M.)] is a nonzero 
ideal in sl(n, R) and sZ(n,R) is a simple Lie algebra. Since dimsZ(n,R) = ?i 2 — 1 and 
B k sl(n,M), dim I > n 2 . Thus t = gl(n,M). □ 



The next lemma is used in the proof of Theorem 13. II in the next section. Let C be 
a Lie algebra, A E C, and let JC be a subspace of L. Recall that the operation cmIa is 
defined as ad.A(B) := [A, B], and the normalizer of K. is 

JV £ (/C) = {4 : [A, B]eK for all B E K). 

It follows from the Jacobi identity that Nc(IC) is a subalgebra of K [321 p. 7]. 
Lemma 2.3. Let A,L € % n (R). Assume (iA,iL)i.j — u(n) and define 

S := B pan({ad*i • ■ ■ ad^" 1 adfe [U, ii]}), (4) 

where s and ki,...,k s are nonnegative integers. Then S = su(n). 

Proof. First note that [iA, iL] ^ because we have assumed that iA and iL generate 
u(n). Clearly iA,iL E N u ( n )(S). Since A r u („)(5) is a subalgebra of u(n) and iA 
and iL generate u(n), N u r n \(S) = u(n). Thus S is an ideal of u(n). Notice that 
iS C su(n) since [iA,iL] is skew-Hcrmitian with zero trace and iA and iL are skcw- 
Hcrmitian. Since S is an ideal of u(n), S is an ideal of su(n), and S ^ {0}. Since 
su(n) is a simple Lie algebra, by definition it has only the trivial ideals {0} and su(n). 
Therefore S = su(n). □ 

For A E W n (R) and Z = {zi, . . . , z s } C R™, the real Lie algebra generated by A 
and Z is defined as 

C(A,Z) := (A,z lZl T ,...,z s z s T ) M . (5) 

3 Controllability and walk matrices 

For A E ■%„(R) and Z = {zi, . . . , z s } C R", the extended walk matrix of A and Z is 
the n x (ns) real matrix 

W(A,Z) := [z^Azi,...,^"-^!,...^,,^,...,^"-^,]. (6) 

A special case is when Z = Zs := {ej : j E S} (with denoting the j-th standard 
basis vector) for some subset S C V(G) for a graph G and A is the adjacency matrix 
Aq of the graph. In this case, the relevant walk matrix is W(Aq, Zs). 

For s — 1 the connection between the walk matrix W(A, Z) in ^ and the 
Lie algebra C(A,Z) in ([5]) was studied in [10]. It was shown [TOl Lemma 1] that 
rank VF(^4, {z}) = n implies C(A, {z}) = gZ(n, R), or equivalently, {iA, izz T )r.j = u(n) 
(cf. Proposition 12. ip . The next theorem states that the converse of this result is also 
true. 

Theorem 3.1. Consider a matrix A in "H„(R) and a vector z E R™ xn . Then, 
(iA, izz T )r. ) .i = u{n) (or equivalently C{ A, {z}) = gl(n,M)) implies that r&nkW (A, {z}) = 
n. 
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Proof. The equivalence of the hypotheses is justified by Proposition l2.fi The result is 
clear if n = f , so assume n > 2. We use a contradiction argument. Assume the rank 
of the walk matrix W(A, {z}) is less than n but (iA, iL)[.,.] = u(n), where L := zz T . 
There exists a vector xgC" such that x*W(A, {z}) = 0. Consider the rank f matrix 
D := xx*. We claim that D commutes with every matrix in S, where S is as in (j4|) . To 
see this, notice that from (|4"|). all elements in S are linear combinations of monomials 
of the form M = A kl L k2 A k3 ■ ■ ■ L kp - 1 A kp , for some p > 1, k 3 > 0, and L appearing at 
least once with exponent greater than zero. When multiplying D with M, with D on 
the left, write M as A kl LY for some matrix Y, so we have 

DM = DA kl LY = xx.* A kl zz*Y = 0, (7) 

which follows immediately from the condition x*W(A, {z}) = for n — 1 > k\ > 0, 
and by using the Cayley-Hamilton theorem for k\ > n. Analogously, when multiplying 
D on the right of M, we write M as QLA kp , for some matrix Q, and we have 

MD = QLA kp D = Qzz*A kp xx* = 0, (8) 

since x*A kp z = also implies z*A kp x = 0. Therefore D commutes with all elements 
of S. 

Observe that since su(n) is simple, su(n) is an irreducible representation of su(n). 
Therefore, since D commutes with all elements of 5, it follows from Schur's Lemma 
that D must be a scalar multiple of the identity [HI p. 26]. However this is not 
possible since D has rank 1. This gives the desired contradiction and thus completes 
the proof. □ 

We study the generalization of this result to multiple vectors (s > 1) but for 
matrices A and vectors Zi,...,z s related to a connected graph G. In particular, 
Q(A) = G, all nonzero off-diagonal entries of A have the same sign, and , . . . , ej s 
will be the characteristic vectors associated to a subset S of the vertices. In the 
next section we will relate this to the zero forcing property of the set S. In the 
context of graphs, it is important tojeonsider multiple vectors because if G is a graph 
and rankAc < \G\ — 2, then rank W{Aq, {z}) < n for any one vector z. On the 
other hand we will see that if S is a zero forcing set for G and G(A) = G, then 
C(A, {e 3 :jeS}) = TL n (R) (see Theorem gj] below) . 

The next definition extends the definition given in [3] (and implicitly in [10] ) of an 
associative algebra that links the walk matrix and controllability. 

Definition 3.2. For A 6 H n (R) and Z = {z 1: . . . ,z s } C K™, define 

P(A,Z) := {A m z, k <L j T A 1 : 1 < k,j < s,0 < m,£<n-l}. 

Remark 3.3. For A s 7^„(R) and Z = {zi,...,z s } C R™, the associative algebra 
generated by P(A, Z) is equal to span P(A, Z), because 

(A m z k Zj T A e ){A a z p z q T A h ) = (z/A e+9 z p )A m z k z q T A h and z/^+% G R. 

Lemma 3.4. For A e H„(R) and Z = {z x , . . . ,z s } C R", r&iikW{A, Z) = n if and 
only if span P (A, Z) = R nx ™. @ 

4 As a vector space, R nxrl is the same as gl(n,M.). We use the latter notation when we want to stress 
the Lie algebra structure on gl(n, R). 
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Proof. Clearly ra,nkW(A,Z) = n if and only if range W(A, Z) = K™. First assume 
ia.nkW(A, Z) = n. For any matrix M G M. nxn with rank M = r, there exist vectors 
x(«),y(«), (? = l,...r, such that M = Eg=i x(?) y (5f)T Since range W(A,Z) = K™, 
each x^ is expressible as a linear combination of the columns of W(A,Z), i.e., as 
a linear combination of vectors of the form A m Zk, and similarly for y^ q \ Thus each 
x WyW : anc j hence M, is expressible as a linear combination of A m z k Zj T A e . Thus 
the matrices of the form A m z k Zj T A 1 span R nx ". 

For the converse, observe that if B = {b^ . . . , b r } is a basis for range VF(^4, Z), 
then 

spanP(A, Z) = spandb^b., 71 : 1 < k,j < r}). 

If n > r = rankW(A, Z), then dim span P (A, Z) < r 2 < n 2 = dimR nxn , so the 
matrices in P(A, Z) cannot span R nx ™. □ 

The distance between two distinct vertices u and v of a connected graph G, denoted 
by d(u, v) is the minimum number of edges in a path from u to v. 

Lemma 3.5. Let A G T-L n (M.) such that G(A) is connected and all nonzero off-diagonal 
entries of A have the same sign. IfkJ E {1, . . . ,n} andk^ j, then (A d ^) kj 0. 

Proof. Let d := d(k,j). The entry (A d ) kj is a sum of terms which are each the 
product of d nonzero entries of A. Since d is the distance between k and j, only 
off-diagonal entries can appear in this product. Thus every term has the same sign 
and {A d ) k] ^0. □ 



Lemma 3.6. Let A £ H n {M.) be such that Q(A) is connected and all nonzero off- 
diagonal entries of A have the same sign. Let S C {1, . . . , n) and Z = {ej : j G S} be 
the subset of standard basis vectors. Then spanP(A, Z) C C(A, Z). 

Proof. The proof of Lemma 1 in |10| shows that for any real symmetric matrix A and 
vector z, A m zz T A i G C(A, {z}) for all m, £ G {0, . . . , n — 1}. Applying this, we obtain 
that A m e j e J T A e G C(A, Z) for all j G {1, . . . , s}, m, I G {0, . . . , n - 1}. The result will 
follow if we are able to show that A m ek^j T A 1 G C(A, Z) for all k,j G {1, . . . , s}, m, £ G 
{0, . . . , n — 1}, with k different from j. 

Consider the distance d(k,j) between the nodes k and j in G(A), which is < n — 1 
because Q{A) is connected. From the fact that both efcejT and A^'^ejeJ are in 
C(A, Z), we have in C(A, Z), 

[e k e k T ,A d ^)e jej T ] =e k e^A^e j e j T -A d ^e j e j T e h e k T = (e k T A^e^ef . 

It follows from Lemma 1531 that e k T A d ^ej ^ 0, and so e k ej T G C(A, Z). 
Then 

[^ m e fc e fc T ,e fe ej T ] = A m e k e k T e k e j T - e k ej T A m e k e k T 
= A m e k e, T - (e/A m e k )e k e k T . 

So, A m e k ej T G C(A,Z). Similarly, e k e 3 T A^ G C(A,Z). Finally, 

[A m e k e k T ,e kej T A e ] = A m e k e k T e k e 3 T A e - e k e 3 T A m+e e k e k T 
= A m e fc e/^-(e/A m +^)e fceib T . 

So, A m e k e j T A i e £{A,Z). □ 
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The following theorem establishes the connection between quantum Lie algebraic 
controllability and the rank condition for an extended walk matrix modeled on a 
graph. 

Theorem 3.7. Let A G H n (M) such that Q(A) is connected and all the nonzero off- 
diagonal elements of A have the same sign. Let S C {1, . . . , n} and Z = {e^ : j G S} 
be a subset of standard basis vectors. Then rank VT(j4, Z) = n if and only if £(A, Z) = 
gl(n,R). 

Proof. By Lemma[320 spanP(A, Z) = R nx " if and only if iankW(A,Z) = n, so it 
suffices to show that spanP(A, Z) = R nxn if and only if C(A,Z) = gl(n,R). By 
Lemma EH span P(A,Z) C C(A,Z), so span P{A,Z) = K ?IX " implies C(A, Z) = 
gl(n,R). For the converse, assume £(A, Z) = gl(n, R). Then, by Lemma l2~2"1 £ = 
gl(n,M.), where £ is the smallest ideal of £(A, Z) that contains ejej T , j = 1, . . . , s. It 
is clear that £ C spanP(A, Z), so spanP(A, Z) = R nxn . □ 



Corollary 3.8. Let A G H n (M.) such that 9(A) is connected and all the nonzero off- 
diagonal elements of A have the same sign, and let S C {1, . . . , n}. Then rank {e^ 
j G S}) = n if and only if (iA, {iejej T : j G >?}}[.,.] = u(n), i.e., the quantum system 



associated with the Hamiltonians iA and iejej 



3 = 1, 



is controllable. 



Observe that for any connected graph G, the adjacency matrix Aq and the Lapla- 
cian matrix Lq satisfy the hypotheses of Theorem 13 . 71 and Corollary 13. 81 

The result of (TO] for the case s = 1 showing that rank W(A, {z}) = n implies 
(iA,izz T )[.^ = u(n), (and the converse proved in Theorem 13.11 in this paper) were 
proved in reference to systems on graphs. The proofs however go through for an 
arbitrary symmetric matrix A and vector z. It is natural to ask whether the conditions 
on the matrix A that we have used in Theorem 13.71 are really necessary. To this 
purpose, we can observe that the result is not true if we give up either of the hypotheses 
that 1) G(A) is connected or 2) the off-diagonal entries of A have the same sign, as 
shown in the next two examples. 

Example 3.9. To see the necessity of assuming that Q(A) is connected, consider a 

'At " 



block diagonal matrix A — 

U Ai 

sions ni and n.2, respectively, with n\ 



with Ai and A 2 symmetric matrices of dimen- 



Ti2 = n, and z\ and Z2 two vectors that have 
zeros in the last n 2 or first ni entries, respectively, and such that the corresponding 
matrices VT(Ai,{zi}) and VF(A2,{z2}) have ranks n\ and ri2, respectively. In this 
case, the walk matrix W(A, {zi, Z2}) has rank n, but the Lie algebra generated by A, 
zxzj, and Z2z|" contains only block diagonal matrices. 

Example 3.10. To see the necessity of assuming that all nonzero off-diagonal entries 

TO 1 ll 



of A have the same sign, consider A 



and Z = {e!,e 3 }. It is 



10-10 
0-101 
10 10 

straightforward to verify that the walk matrix {e!,e 3 }) has rank 4. However, 

rank£(A, {e!,e 3 }) < 8, as can be seen as follows. Let 
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"1 o ol [0 0] [0 O Ol [0 0" 
, 0000 0100 0000 1000 
: - s P ai K OOOO'OOOO'OOIO'OOOO' 
_o o o oj |_0 lj |_o o o oj |_1 

"0 1 1] [0 Ol [0 Ol [0 0' 

0000 0000 00 1 0001, 

0000' 010-1 '00 ' o o 

o o o oj [o o J [o -1 oj |_0 1 

Since [B, C] G C for all B, C € C, C is a Lie subalgebra of gl(4, R). Clearly dim C < 8 
and C(A, {ei, e 3 }) C C. 



4 Zero forcing and controllability 

The neighborhood of v G V(G) is N(v) = {w G V(G) : w is adjacent to v}. 

Theorem 4.1. Let A £ "H„(R) suc/i </iai 0(A) is connected and all the nonzero off- 
diagonal entries of A have the same sign. Let V := {1, 2, . . . , n} be the set of vertices 
for Q{A), and S C.V be a zero forcing set of Q{A). Then 

M <>,<'/' :./ £ S}) = gl(n,R). 

Proof. After a (possibly empty) sequence of forces, denote by T the set of currently 
black vertices, and assume that for all k G T, ekek T G C := C(A,{ejej T : j G 5}). 
The hypotheses of Lemma \3. 61 are satisfied for Z = {e^ : j G T}, so for all k,£ G T, 
e k e e T G C. 

If T V, then there is a vertex u G T that has a unique neighbor w outside T. 
For that u we have 

[^u^tz ; A\ ^ y o urn ( y G u e rn e m e u ), 

mgW(u) 

where a um := a um if u < m and a um :— a mu if m < u. For all m G N(u) such that 
m ^ w, m G T, so e u e m T ,e m e u T G C. Thus e„e^ — e^e^ G C. Since 

r T T 1*1 T T 

|e u e u , s u b w g w g u J e^e.^ -I- e^e^ , 
e w e^,e w e^ G C. Then 

r T Ti T T 

so e^e^ G C. Since 5 is a zero forcing set, we obtain e£e m T G C for all £,m G V = 
{1, . . . , n}, and thus we conclude that £ = gl(n, K). □ 

Applying Proposition 12 . 1 1 we obtain the next corollary. 

Corollary 4.2. 7/G is a connected graph, A G "H n (IR.), aZZ £/ie nonzero off-diagonal en- 
tries of A have the same sign, and S C V is a zero forcing set of G, then (iA, {iejej T : 
j G = u(n) and the corresponding quantum system is controllable. 

Note that the converse of Theorem 14. II is false. 
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Example 4.3. Consider the path on four vertices P4 with the vertices numbered in 
order. The set {e2} is not a zero forcing set for P4. However, 



W(A Pi ,{e 2 }) = 



10 2 
10 2 
10 3 
10 



and rankW(A Pi ,{e 2 }) = 4, 



so C(A Pi ,{e 2 }) = gl{n,R) by Theorem GEO 



5 Conclusion 

Motivated by the control and dynamics of systems modeled on networks both classical 
and quantum, we have established a connection between various tests of controllability 
and the notion of zero forcing in graph theory. Lie algebraic quantum controllability 
is necessary and sufficient for linear (Kalman-like) controllability of an associated 
system and both notions are implied by the zero forcing property of the associated set 
of vertices. Linear systems have a very well developed theory [T5] and it is an open 
question to investigate to what extent this analogy can be further used to discover 
properties of quantum systems and systems on networks. 
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